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Abstract In two space dimensions, the parabolic-parabolic Keller-Segel system shares 
many properties with the paraboHc-elhptic Keller-Segel system. In particular, solutions 
globally exist in both cases as long as their mass is less than 8 vr. However, this threshold 
is not as clear in the parabolic-parabolic case as it is in the parabolic-elliptic case, in 
which solutions with mass above 8 tt always blow up. Here we study forward self-similar 
solutions of the parabolic-parabolic Keller-Segel system and prove that, in some cases, 
such solutions globally exist even if their total mass is above 8 7r, which is forbidden in 
the parabolic-elliptic case. 
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1 Introduction 

The Keller-Segel model has been widely studied for almost forty years. It models the 
behavior of a slime mold of myxamoebae, Dictyostelium Discoideum, which have the 
peculiarity of organizing themselves to form aggregates by moving towards regions of 
a higher concentration of a chemoattractant. This chemoattractant, the cyclic adeno- 
sine monophosphate, is secreted by the amoebae themselves when they are lacking 
of nutrients. The Keller-Segel model is considered as a prototypical model for pat- 
tern formation in chemotaxis, and has attracted a lot of attention as a test case for 
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more complex taxis phenomena driven by chemical substances. See [9l ll0l[Tl1ll7lll8j for 

further references. 

The simplest version of the model is made of two equations, one for the density 
of the amoebae and another one for the density of the chemoattractant. Both are 
parabolic, although an even simpler version has been widely considered by neglecting 
the time-dependence of the density of the chemoattractant. We shall refer to the com- 
plete version of the model as the parabolic-parabolic model, and to the latter as the 
parabolic-elliptic Keller-Segel model. 

The precise dependence of the diffusion coefficients and of the chemosensitivity 
parameter depends on the context. All variants of the model involve diffusions in the 
equations for the density of the amoebae and for the density of the chemoattractant. 
The coupling is due to the fact that amoebae move according to the gradient of the 
chemoattractant, and that the emission of the chemoattractant is proportional to the 
density of amoebae. A crude insight into the main features of the model can be gained 
from the simplest case, that is when the nonlinear term in the equation is quadratic, 
but more realistic models should probably involve more complex nonlinearities. 

Since the slime mold usually moves over a planar substrate, it makes sense to 
consider two-dimensional geometries. In some cases, boundary effects are important, 
but they are out of the purpose of this paper and we shall therefore assume that the 
model is set on the two-dimensional Euclidean plane. 

In the parabolic-elliptic model, there is a critical mass, 8 tt after a proper adimen- 
sionalization, whose role is now rather well understood; see [8][5]. Below such a mass, the 
diffusion predominates, in the sense that amoebae are unable to emit enough chemoat- 
tractant to aggregate. On large times, the population diffuses and locally vanishes, 
although the behavior significantly differs from a pure diffusion. Above Svr, at least 
one singuiarity appears in finite time, which is interpreted as the occurrence of an 
aggregate. 

Since singularities are local, it is widely believed that 8 tt should also be a threshold 
between the diffusion dominated regime and the regime of aggregation also in the 
parabolic-parabolic model. This is certainly the case in some sense, for appropriate 
initial data, but the situation is not as simple as in the parabolic-elliptic case. It turns 
out that if, initially, the population of amoebae is scattered enough, and for a well 
chosen initial distribution of the chemoattractant, there are solutions for which the 
diffusion predominates for large times, even for masses larger than 8 tt. It is the purpose 
of this paper to establish such a fact, for a special class of solutions and in a certain 
range of the parameters of the model. 

In this paper, we consider the parabolic-parabolic Keller-Segel model 



for the densities n and c of, respectively, microorganisms {e.g. amoebae) and diffusing 
chemicals that they are secreting. Interesting mathematical questions are related to 
qualitative properties of problem ([lj-((2ll such as global in time existence versus finite 
time blowup of solutions describing chemotactic concentration phenomena. After the 
pioneering works of Keller and Segel, a huge literature has dealt with the mathematical 
modelling of chemotaxis and its analysis. We recommend the reading of [9] for a recent 
review from both biological and mathematical points of view. 



nt — An — V ■ (n Vc) , 
T ct ~ Ac + n , 



(1) 
(2) 
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We shall consider the Keller-Segel system (IT|-(l2} for any t > 0, x £ M.^ , supple- 
mented with initial conditions uq and cg • From now on we shall assume that uq and cq 
are nonnegative and that ng is integrable on R^. As a consequence, for solutions with 
sufficiently fast decay at infinity, the total mass is conserved, i.e., 

M := n{t,x)dx— / 710(2;) da; 

does not depend on t. 

Throughout the paper, r is a nonnegative parameter taking into account the differ- 
ence of the time scales of the diffusive processes undergone by 71 and c. The qualitative 
properties of n and c (such as the asymptotic behavior for large values of t) depend 
on r and the stability of system (Il|-((2} with respect to r is expected, i.e. solutions 
of the parabolic-parabolic Keller-Segel system are expected to converge to those of 
parabolic-elliptic system when r \ 0. This has been recently proved, at least for solu- 
tions with a suitably small mass M, in p2]. Here, we are interested in the differences 
between the parabolic-elliptic Keller-Segel system (r — 0) and the parabolic-parabolic 
Keller-Segel system (r > 0). Known results are briefly summarized as follows. 

When r = 0in(l2]),M = 87risa threshold for existence versus blowup of the 
solution of (Il]|-((2]), see [8]|5l[7]. Solutions globally exist for M <8n, while explosion in 
finite time may occur if M > 8 tt. In the critical case M = 8 tt, the solutions are known 
to be global in time but the density grows and mass concentration occurs in infinite 
time; see [3l|4]. 

For r > 0, according to JT, solutions globally exist for any M < 8n. However, it has 
not yet been proved that explosion occurs in finite time as soon as M > 8 tt, for instance 
under some additional assumptions like a smallness condition on |a;|^ 710(2;) da;. If 
M = 8tt, there is an infinite number of steady states (see [3]), but no other result is 
available, apart from self-similar solutions. 

Motivated by this lack of results for (Il])-([2} , this paper deals with the existence of 
positive forward self-similar solutions of ([T])-((2|, i.e., solutions which can be written as 



n(t, x) — — u ( and c(t,x) = v(-^] 

t \VtJ \VtJ 



(3) 



with a large total mass (that is, larger than Svr). Indeed, since we are dealing with 
the two-dimensional case, any self-similar solution n in L^(R^) preserves mass, i.e., for 
each t > 



/ n{t,x)dx= / M(C)dC = M. 



Therefore, for any given r > 0, we are interested in the optimal range of M for the 
existence of such solutions, and in uniqueness or multiplicity issues for a given M in 
the optimal range. Actually, our goal is double. The main one is to prove the above 
mentioned existence result. Second, we will give an as complete as possible review of 
the numerous existing results on the topic and also simplified, new proofs of them. For 
this reason, the remainder of the introduction will be primarily devoted to the state of 
the art on self-similar solutions. 



Self-similar solutions can be obtained through various approaches. The first method 
for the study of self-similar solutions (see for example [T] and the references therein) 
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amounts to look for mild solutions of I!}-©, that is, solutions of 

n{t, .) = e(*-*°)^ n{to, •) - j^' (Ve^*-^)^) • {n{s, •) \/c{s, •)) ds , 

c{t, ■) = e ^ c(to, ■) + -"/ n{s, •) ds , 

Jto 

for any t > to > 0- Roughly speaking, such self-similar solutions are obtained by a fixed 
point theorem. However, smallness conditions on the initial data are required in order 
to apply a contraction mapping principle; see [13] , where this method has been applied 
to I!])-® with r = 1. Therefore, covering the whole range of masses for which solutions 
exists seems out of reach in this setting. 

Alternatively, one can prove the existence of self-similar solutions through the direct 
analysis of the elliptic system satisfied by {u,v), i.e., 

Au-V ■ ^uVv-^^v^ = 0, (4) 

Av+^^-\/v + u = 0, (5) 

where ^ — xjyft and the differential operators in Q-® are taken with respect to ^. In 
this case, a natural functional space to be considered for both u and v is the subspace 
Cq(R^) of functions in the space C^(R^) such that 

lim uii) = and lim v(f) = . 
For such classical solutions, equation Q can be written equivalently as either 

u V (log u-v+^-^ ] = , 



V- e"e 



.^.-|«lV4. 



Then, using the fact that u, v, and consequently IVnl are bounded, it has been proved 
in that there exists a constant a such that 

^^(0 ='Te^(«)e-^ (6) 

for any ^ £ R^. Since u is positive by the maximum principle, it follows that a is 
positive. As a consequence, u £ L^(R^), and the stationary system ©-([5]) reduces to 
a family of nonlinear elliptic equations for v, namely 

T i;p 
Av + -^-Vv + ae''ei-=Q, (7) 

parametrized by a > 0. Again by the maximum principle applied to ([71, the following 
upper bound for v can be proved 

i;(e)<Ce-™"{l'">^, (8) 

where C is any positive constant such that C min{l, t} > a e"""°° ; see for instance [15j . 
Therefore, -u G L^(R^) holds true for any solution of in Co(R^). 
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The range of M for which self-similar solutions exist in Cq (R^ ) gives an indication 
on the range of M for which some solutions of H}-© may globally exist. Self-similar 
solutions indeed provide explicit examples of global solutions, even with smooth initial 
data, up to a time-shift: take for instance u and v as the initial data for H}-©. 
Moreover, if self-similar solutions describe the asymptotic behavior of any solution of 
(IT])-([2]) under appropriate conditions on initial data, then the ranges of global existence 
of solutions should be exactly the same. This property has been established in [5] for 
T = 0. In the case r > 0, this might not be as simple as in the case r = if one can 
prove that blowup may occur for any M > 8n. However, at least for initial data close 
enough to u and v, one can expect that the ranges of global existence are the same. 

In view of our main goal, we are actually more interested in parametrizing the set 
of Co(R^) self-similar solutions in terms of mass rather than in terms of a. This is 
possible using in (O the relation 

M = a e"(«) e — ^ d^ . (9) 

However, by doing that, equation ([7]) becomes nonlocal, as was the original system 
([T])-((2|, and the problem is definitely more difficult to handle. Another not less im- 
portant reason to consider a different but equivalent formulation of problem ©-((Si is 
that the correspondence between a and Af is not clear due to the lack of uniqueness 
of solutions to Q; see Remark [T] at the end of Section [2] 

For the sake of completeness, we have to say that equation (0, written as 

V • (e^^^^^ Vv^ + cr e" e^l^l"" = , 

has been studied using variational methods in [131119) . The weighted functional space 
//^^(R'^; exp(^I^P) d^) is then natural, but working in this space introduces a condition 
on the values of r, which have to be in the interval (0, 2). Under such a restriction, 
it has been established that solutions exist if < a < a*, for some a* > 0. These 
solutions are positive and belong to Cq(R^), but due to the restriction on r, one has 
to look for alternative approaches. 

Another important and useful result has been obtained in [15 using the moving 
planes technique: any positive solution v G Cg(R^) of ((TJ must be radially symmetric. 
As a consequence, system @-(0 reduces to the ODE system 

u— uv' + —ru = 0, (10) 

+ / + u = 0, (11) 

where u and v are considered as functions of the radial variable J" = |CI oiAy. Equations 
©-(O then become 

/ \ v(r) — 

u[rj = a e ^ ' e , 

«"+ (^i + Ir^^;'+ae''e-'■'/4 = 0. (12) 

Equation (|12|) has been studied in |12II15) . More specifically, the authors proved 
in [12j the existence of a positive decreasing solution of (1121) endowed with the initial 
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and integrability conditions 

/■OO 

v' (0) = and / rv{r)dr<oo, (13) 
Jo 

for any r > and a > such that cr -^rx < 1/e (see Remark [2]). However, such a 
condition does not determine the optimal range neither for the parameter a nor for M. 
It is worth noticing that the boundary conditions (|13p and the following ones, 

«'(0) = and limu(r) = 0, (14) 

r — >oo 

are equivalent for classical decreasing solutions. Indeed, (|13[l implies (|14|l and the con- 
verse holds true by ([8]). Using (|14|) . equation (|12p turns out to be equivalent to 

™" + (^i + ^ ^i;' + e"" e-'-'/^ = , (15) 

iu'(0)=0 and ™(0) = s , (16) 

for some shooting parameter s G R. Indeed, if w(r; s) is a classical solution of (|15|l - 
(|16|1 for a given s G R, then w(cx); s) = lim^— »oo M(r; s) exists and is finite and v{r) = 
w(r; s) — ui(oo;s) is a classical solution of p2p - (ll4|l with a = e"'^""'*''. Conversely, if 
u is a classical solution of p2|) - (ll4|l . then w{r;s) = v{r) + log cr is a classical solution 
of ([T5l)"(IT6| with s = v{0) + log cr and again a = e'^^°°'^'> holds true. It follows that 
all solutions of (|12|| - (|14p can be parametrized in terms of s. See [T^ for more details. 
Using this equivalence, the authors of [15] analyze the structure of the set of solutions 
of p2p - (|14l) seen as a one-parameter family; see Remark [T] at the end of Section (2] 
Computations presented in Figs. 1 have been based on this parametrization of the 
solution set. 

Last but not least, the parametrization of the solutions of (|15 [l -(|16 p in terms of s 
allows us to parametrize the total mass Af in term of s by 

roc 2 / 

M(s) = 2 7r/ e'^^'^-^'K-'' /*rdr. (17) 
Jo 

But again, this does not provide an explicit computation for the optimal range of M. 
Computations presented in Fig. 2 (left) have also been based on this parametrization 
of M. 

Being this the state of the art, we will establish that the formulation of system 
pop - (|lip in terms of cumulated densities is better adapted to the qualitative descrip- 
tion of u and v. This is a classical technique used previously, for example, in the context 
of the parabolic-elliptic Keller-Segel system and astrophysical models; see [T]|3] and fur- 
ther references therein. For r > 0, many qualitative properties of the solutions can still 
be proved in this framework. These will allow us to build positive forward self-similar 
solutions of HI-® satisfying which have an arbitrarily large mass when r is large 
enough. The obtained results are summarized in the theorem below. One may inter- 
pret it by saying that the diffusion of c described by ((2} for positive large r and some 
M > 8 TT may prevent the blowup of the solutions of the parabolic-parabolic Keller- 
Segel system. This is a major difference with the parabolic-elliptic case r = 0, for which 
the response of c to the variations of n being instantaneous, any smooth solution with 
mass M > 8 TT must concentrate and blow up in finite time. 
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Theorem 1 For any M > 0, there exists some f{M) > such that for any r > f{M) 
there is at least one solution {u,v) of ((4|-((5]) in (Co(K^))^ with u > of mass M 
and V > 0. If M < 8tt, t{M) = 0. If M > 8tt, f{M) is positive and there are 
at least two solutions, except for the maximal possible value of M . All solutions are 
radial, nonincreasing , with fast decay at infinity, and hence attain their maximum at 
X = 0. They are uniquely determined by a := m(0)/2, which in turn uniquely determines 
M = M{a,T). Moreover, lima^oo M{a,T) = 8tt, while, as a —> oo, the corresponding 
solution u concentrates into a Dirac delta distribution, up to the factor Stt, and w(0) = 
II^IIl°°(R2) becomes arbitrarily large. 

This paper is organized as follows. We shall first establish the main a priori es- 
timates for Theorem [1] in the next section. The framework of cumulated densities is 
developed in Section [S] which also contains more detailed statements than the ones of 
Theorem [1] The remaining a priori estimates and proofs are given in Sections |4] and [S] 
respectively. Section|6]is devoted to some numerical results and Section[7]to concluding 
remarks. 



2 Large mass positive forward self-similar solutions 

Before restating the question of self-similar solutions in terms of cumulated densities, let 
us establish the key a priori estimate for Theorem [TJ which proves that these solutions 
may have an arbitrary large mass when r is large enough. This result is entirely new. 
Such an estimate can be obtained both from equation (I12|l and from the cumulated 
densities formulation. In this section, we shall establish this a priori estimate in the 
first setting. It will be translated in the cumulated densities framework in Section 

From now on, we shall parametrize M in term of a and r, i.e. M = M{a,T), 
where a = u{0)/2 will be the shooting parameter in the cumulated densities shooting 
problem, see and ([33l)-([34l) below. 

A positive classical solution v of (|12|) . (|14|) solves 

(re-'^'/%')'-fare(--l)'^'/4^-^0, 
which, after an integration on (0,r), gives 

1/ N _Tr^/4 (t-1)2^/4 v(z) , 

V (r) — e ' / e^ ' ' e ^ ' z dz . (18) 

r Jo 

As a consequence, v' is nonpositive, so that v{z) < v{0) for any z > and, for r 7^ 1, 

(19) 




Hence, after one more integration of (|19|l on (0, 00), we get, for any r 7^ 1, 

«(0) < ae''^°^ I{t) with /(r) := . (20) 

r — 1 



We observe that 
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Actually, it is easy to check that estimate (|20p holds true also for r = 1 with J(l) = 1. 
Since from (|6]) we have 

o-e^f") =m(0) = 2o, (21) 
it has been proved that for each r > 0, 

= lim v(z) < v(r) < v(0) < 2a/(r) (22) 

Z — ^CJO 

for any r G R+. On the other hand, by ((9)1, (I14|) and (1221) . mass can be estimated for 
any positive a and r by 

oo 2 /"OO 2 

M = 2 TT (T / e^'-''-' e"*^ r dr > 2 n cr e^^^rdr = 4:na>8na e"^" ^^^^ 



_ _ (23) 

using (|2ip . As a function of a, M{a,r) := 8 7rae 2a/(r) ^^jjigygg j^g maximum at 
a*(r) := jT^i which proves that M = M{a,T) verifies for each r > 

4 7r 

max M(a,T) > Af(a*(r),r) 



a>0 ^ ' ' - '-^""-vw, ■ / gj(^) , 

and it is clear that the right hand side can be made arbitrarily large for r large enough. 
Hence, the corresponding density u{r) — ae^^'^'^e~^ 1^ has mass M > Stt if j-j^ ^ 

8 TT, that is for any r > r with f such that /(r) = j^, i.e. f x 16.1109. Also observe that 
for any r > f the density u corresponding to a = a*(T) satisfies u{0) = 2a*(r) > 2e. 
Finally, using 11(2:) > v{r) in (|18p and integrating the inequality on (0, 00), one obtains 



-lim^-.-^ 
limr^oo e""'-'^' = 1, we obtain that 



006 < —aI{T), for any r > 0. As a consequence, using (|2ip and 



l_e"(0) <-a/(r)e''(°) =-2a/(r). 

This gives the estimate 

«(0) > log(2a/(r) + l), (24) 

which implies that v{0) becomes arbitrarily large as a — + cxj, for any r > 0. 

Estimates (I23p and (|24[l can be read also as lower and upper bounds for a = 
2ae""(°), namely 

o -2a/(r) ^ ^ ■ ( M 2 a \ 

2ae ^ ' < a < mm < — , — > , (25) 

- - \47r ' 2a/(r) + l J 

hence showing that a takes arbitrarily large values for r large enough. 



Remark 1 Estimates (|25|l on cr are new. The authors of [15) analyzed the map s ^ 
(t(s), where s is the shooting parameter defined in (|16p . and they proved that it is 
a continuous map from R into IR+ with lims^±oo o"(s) = 0. Therefore, a must be 
bounded for any fixed r by a* = cr(s*), for some s* G R, and problem (|12p - (ll4p 
admits no solution for cr > a*, at least one solution for a = g* and finally (at least) 
two distinct solutions for < a < a*. However, estimates on o (or a*) were missing. 

Remark 2 Estimate (|20|l says that, for any fixed cr > and r > 0, ^(O) satisfies 

w(0) - aJ(r)e"(") < 0. 

Since the function x 1-^ x — a I{T)e^ is strictly concave and attains the maximum in 
X — — log(c7 J(r)), we deduce that whenever u I{t) < 1/e, there exists an open interval 
J C K+ of non existence of solutions of (|12p satisfying (|14p . with v{0) G J. On the 
other hand if o- I{t) > 1/e, the above inequality induces no restriction on t;(0). 
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3 Cumulated densities and main results 

Let us introduce the cumulated densities formulation of the parabolic-parabolic Keller- 
Segel model as in [1], in terms of the functions u and v which solve problem (I10p - (|lip . 
by defining 

1 f 

:= 75— / ru{r)dr, 
1 f 

V'C?/) := 75— / <^i= rv{r)Ar. 



Using the relations 



<f>'{y) = ^ " i^/y) and (t)"(y) = -^-j= u {yfjj) , (26) 

{y) = \'" iVv) and = j^"' , 

it follows from (|10p - (|lip that the cumulated densities cf) and tp solve the second order 
ODE system 

0" + i0'-2 0V' = O, (27) 
Ayxp" + Ty'4)' -Tit) + (I) = Q, (28) 

where has been multiplied by r and integrated on (0, y^) . Observing that equa- 
tion (1281) can be written as 



and defining S{y) := 4(i/;(y) - y^'iy))' = -4yi/'"(i/) = -y^v'{y/y) as in [IKTS], 
system (|27p - (|28|) becomes, after a differentiation of (|28p with respect to y, a first order 
system in the {(j)' , S) variables 

cb" + ^^' + ^c^'S = 0, (29) 
4 2y 

S' + ^S = ^'. (30) 

The last formulation of the ODE system can be equivalently written as a single integro- 
differential equation, hence nonlocal, for 0', 

<P" + U' + ,!-^'e'^y/^ r e^''W{z)Az = Q, (31) 
4 2y Jo 

since, by (f5IJ)) . 

5(t/) = e-" e" (l,'(z) Az , (32) 



and as a single, local but nonlinear second order ODE for 5", 
^" + i(. + l)5' + ^5 + J-(55' + l5^ 
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which is obtained by differentiating (|30p . We will use in the sequel all these formulations 
in order to get a priori estimates. 

For any positive self-similar solution {u,v) £ (Cg(R^))'^, the natural initial condi- 
tions for (I59l)-((30| are 

(l>{0) = , (l>'{0) = a > and ^(O) = , (33) 

in view of the definition of </!> and of (|32p . Moreover, for any self-similar solution u € 
L^(R^), the corresponding cumulated density (p satisfies the boundary condition 

cPioo) := hm <^(y) = . (34) 



The problem is now formulated in terms of a shooting parameter problem (|29|| - (|30p . 
P3I) . with a new shooting parameter a which is directly related to the concentration of 
the self-similar density u around the origin, since a — it(0)/2. This has been obtained in 
Section[2]and will be made more precise below. Let us observe that the relation between 
a and the shooting parameter s defined in ((16} is 2 a = e", since s — v{0) -I- log a. Thus, 
a one-to-one relation is established between the initial valued problems (|29|l - (|3Up . p3|l 
and (|15p - (|16[) as soon as an existence and uniqueness result is established for one of 
them. Moreover, we have 

^ Jo ^ 

and the boundary condition lim,-^oo v(r) — is equivalent to 

«(0) = ir^d^. (35) 
^ io 2 

We also have: a = limr— »oo u-ir) ~ 21imy^oo 'P'iy) e^^^. Hence we can reparame- 

trize 11(0) and a in terms of a and r. 

The main statements we are going to prove are summarized in the following theo- 
rems. The a priori estimates will be established in Section |3] The proofs will be given 
in Section [S] We shall say that (</>, S) is a positive solution if both (j) and S are positive 
functions. 

Theorem 2 For any (a, r) G E'^ there exists a unique positive solution (</!>, S) of (|29p - 
PU)) . pS)) such that (j) e C^(0,cx)) n C^[Q,oo) and S G C^[0, 00). Moreover, for any 
fixed T > Q, (j) £ C^[0, 00), the maps a £ ^^ (<^, -S") and a G R+ M(a,r) G K-|- 
are continuous and 

I \ ^ M{a,T) 
g(a,T) < '- < f(a,T), 

2 TT 



where 



mm 



{4,4a} tf re (O, i] , 

{4«,|vr2} if re , (36) 



.{4a,§^2^,4(r + l)} if 



r > 1. 
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and 



g{a,T) = < 



'max{4ae-2'^^,iii} ^f re [0,1]., 
,x{4ae-2''^,^} ^f r>l. 



(37) 



For consistency, it is worth noticing that the inequahty g{a,T) < f{a,T) holds for 
aU r > and a > 0. 



Theorem 3 Given any fixed t > 0, for any positive sequence {aj.} such that a^, cxj 
ask ^ oo, there exists a sequence of positive self- similar solutions (ttfe,Ufe) G (Co(R^))^ 
satisfying (SI)-© and Wfe(O) = 2afc, v^(0) = such that 

Uk 8tt 5q as k ^ oo 

in the sense of weak convergence of measures. Moreover, limj-^o^ J^2 vl^. dx = 8 tt and 
limfc_»oo lkfcllL°°(K2) = oo. 

Theorem |3] has aheady been proved in |15l Th. 2, (iii)] using a classical result by 
Brezis and Merle in [6]. However, here we shall give a simplified and quite direct proof 
using the cumulated densities formulation. 

Theorem 4 For any fixed r > there exists M* = M*(r) > Svr such that problem 
(|29[) - p0p with the boundary conditions 

m = 0, hm ^{y) = ^, 5(0) =0, 

has no positive solution (</>, S) £ C^[0, oo) x [0, oo) if M > M* and has at least one 
positive solution {(f), S) £ C^[0, oo) x C^[0, oo) in the following cases: 

(i) M G (0,M*] if M* > Stt, 

(ii) M G (0,M*) tfM* = Stt. 

Moreover, there exist 1/2 < r* < r** such that M* satisfies: M* = 8n if < t < t* 
and M* > 8tv tf t > t** . 

When M* > Stt, there are at least two positive solutions for any M G (87r,Af*). 
When M* = Stt, it is still an open question to decide if there is a positive solution 
{(j),S) G C^[0,oo) X C"'^[0,oo) such that M — M* or to prove a uniqueness result for 
any M G (0,8 7r). 

Remark 3 The estimate r* > 1/2 will be given in Proposition [T] as well as refined 
estimates on Af(a, r). Theoretical results show that r* G (0.5, 16.1109 . . .), see Th. [2j 
(|36|) - (|37p and Sec. [2] while numerical computations suggests that r* G (0.62, 0.64), see 
Fig. 2 (right). Moreover, it is an interesting open question to decide whether r* = r**, 
as again the numerical results suggest, or not. Exact multiplicities of solutions for 
M > 87r are not known in detail either. Let us observe that for t > r* the function 
M(a,r) depends on a in a nonmonotone manner. This is a significant diff'erence with 
the monotone dependence of self-similar solutions of the parabolic-elliptic Keller-Segel 
system (see [3] Sec. 4]). 
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4 Qualitative properties of (p and S 

In the present section we will derive all a priori estimates on (j) and S which are 
necessary to prove Theorems [2] [3] and [l] Some of them are new while other were 
already known. In any case, we shall give a unified and simplified proof of all of them 
in terms of cumulated densities. 



4.1 Preliminary estimates 

Let {u,v) G (Co(R^))^ be a positive solution of Q-Jl} with u G L^(K^). The corre- 
sponding {(j),S) satisfies (l29]|-([30l), (IMl) with a = w(0)/2. Moreover, for any j/ > 0, it 
immediately holds true that: <^ is a positive, strictly increasing and concave function 
on (0, oo) while < S{y) < (f){y) for any y > since S' < (p' on (0, oo) by ((30]). More 
precisely, an integration by parts in (I32p gives 



S[y) = <l>{y) - 7 / ' <l>{z) Az . (38) 







4 

On the other hand, in (|38}, the increasing monotonicity property of (f> gives us 



S{y) > 0(y) - 7 e-^ ''/^ P{y) r e^ d^ = e"^ ''/^ 0(y) , (39) 
4 Jo 

while the decreasing monotonicity property of 0' in (|32p leads to 

S{y) > e-'y/^ p' {y) ''/^ Az = ^<j>' (y) (l-e"^^/'') (40) 

for each y > 0. From (|39|l and (|40| . we get 

J5(y)> (0(y)-0(y)e-^^/4)'. 
Since ^ S = 2<f) — 2S' , the last inequality gives 

S(y)<i0(2/)(l + e-^^/^) 

for each y > 0, which is a better estimate than S < 4> but still not yet satisfactory for 
large y. 

Let us now estimate 4>. Looking closer at system (|29[l - (|3Up , one observes that the 
quantity e^^^ 4> {y) is positive and decreasing. Hence 

l{a,T) ■- lim e^/^ <})' {z) < e^^^ <^'(y) < <j!>'(0) = a, (41) 

z — >oo 

for any y > 0. Notice that l{a,T) — (j/2, which proves that lims_>-_oo o"(s) ~ with 
the notations of Remark [1] Integrating once more the above inequalities on [0, y] we 
have 

4/(a,r) (l - o^y/*^ < <j>{y) < 4a (^1 - e"^/'') . (42) 
In particular, for each r > 0, M[a, r) is finite, 

Z(a,r)<M^<a, (43) 
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and we see that, whatever r is, the shooting parameter a has to be large enough (a > 1) 
in order to obtain a self-similar solution u with mass M > Svr. 

We can improve estimate (I42p as follows. Since limy_>oo 0'(y) = 0, integrating the 
inequality + < on [y, cxj), we get 

X 1 X ^ M(a,T) 
0(j/) + 4 0(2/)> 3^ , 

and therefore, by integrating once more on [0,y], 

0(y)>^(i-e-/^ 

In conclusion, using the previous estimate for (j), we obtain for each y > 

^ (l - e-y'^) < Hy) < min {4a (l - c-/^) , ^} , (44) 



where equality in the minimum is achieved for y = —4 log (l — sttu ) ^ (0,00]. In 
particular, equalities hold in (|44p . i.e. 4>{y) — ^ (l — e~^/^) , if and only if M = 8 tt a, 
in which case y = 00. But since (/>(y) = — e~^/'^) is not a solution of (I29|l - (|30p . 

estimate (|44|) holds true with strict inequalities as well as M < 8n a. 



Coming back to the function S, using estimate (|4ip and identity (|32p . we have 

S{y) < ae-^'^/4 r e^"-^^ '/'^ dz , 
Jo 

i.e., for each y > and r > 0, 

S{y)<ayh{y;T) (45) 

where 



e 



-a/4 if^ = i 



%;^) = <! 4 /^„-i,/4 „-Tw/4\ (46) 

As a consequence, it holds true that 



e-v/-^ ^e-'y''] ifr/1. 



lim S{y) = lim ^ = , 

y — *oo y — >oo y 

S{y)/y is integrable near y — and, using (|40p . ^'(O) — a. 

The above asymptotic behavior of S at infinity, together with the initial condition 
S{0) — 0, allow us to integrate equation (|30p on [0, 00) to obtain 

Mpll=^(^oo) = L rSiy)dy. (47) 
2vr 4 Jo 

Therefore, any appropriate bound for 5 would give a bound for the total mass A/. How- 
ever, let us observe that if we plug estimates (|45p into (|47p . we found again the upper 
bound in (|43p . Finally, thanks to the integrability of S{y)/y near y = 0, equation (|29p 
written as 

and integrated on [0, y] gives the relation 

<t>'{y) = a e-^/4 exp (^-i £ ^ dz) . (48) 
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4.2 Further estimates 

First, let us improve on the lower bound in (|44|) for <j). As far as we know, all estimates 
of this section are new. Using the fact that S < (f> m (I29|) . for y > we have 

4 2y 

After a multiplication by j/, an integration on [0, y\ leads to 
Dividing by (jy^ e we obtain the differential inequality 



-Je--J +-e-->-^e--y^ ^(.) d. . 

Finally, dropping the positive term on the right hand side, and integrating once again 
on [0, y] gives us a lower bound for any r > and a > 0, namely. 



(l + i)e^^/4-l 



for each y > 0. This is, of course, a better estimate than (|44|l but only for y near the 
origin since the inequality S{y) < (jiiv) is a good approximation for y near the origin 
but not for large y. Ifowever, we can now replace (|44p with 

(l - e-/^) . (TTT)^} < ^(y) < {4« (l - e-/^) , ^} . 

(49) 

The maximum on the left hand side of (|49p is achieved by both terms at some y > 
and 

' M(a,r) ( y/4\ V ] _ V 

/ ' (l+i)-"/^-!] (l + i)e!^/4-l 
for each y G [0, y\. Moreover, for any y > y* , we have 

M(a,r) . ^ . *^ y* 4 

if y* is the point where the maximum of y i-^ ^ — j is achieved. 
Next, let us apply estimates (|45ll - (|46)) to (|i8)) . For r / 1, we have 

ry A fy ^ <\ t i py f"^ t 

/ h(z;T)dz=—^ - -^(e-3^)dtdz= — !— / / e-3^dtdz 
Jo T-1 Jo z ' r-lJo Ji 



r-l./itV J r-1 

s{z) 



■dz < 4aI(T) 







with 7(r) = and 



<^'(2/) > ae-^/V^"-^^^). (50) 
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Integrating (I50p on [0, oo), we get the same estimate as in (|23p giving arbitrarily large 
mass M for r large enough, i.e. 

M{a,T) -2a/(T) 

— ^ > 4ae ^ ' . (51) 

For T = 1, since /i(i/; 1) = e"^/*^, one obtains, for aU a > 0, 

M(a,l) 



27r 



> 2 I 1 - e 



Remark 4 The lower bound (|5ip is compatible with the upper bounds for M known 
from [2], I.e. 

:^<4ifrG(0,l/2] , |1 < l^^frG (1/2,1] 

2 TT Z TT O 

M 



and — < min (Itt^t, 4(r + 1)) if r > 1 . 

2 TT L J 



Finally, following 15 , define the new function 

rv 



W{y)~ 0'(z)e^^/'^d2 = e^f/'*S(2/), 





where the second equality follows from (I32p . After a multiplication of (|3ip by e^^^^, 
it is easy to see that W satisfies the initial value problem 

w" + ^w' + — (w^)'e'^y^^ = o, 

4 Ay \ J 

W{0)^0, W'{0) = a. 

Next, a multiplication by y and an integration on [0, y] gives us 



Jo 4 16 Jo 



Dividing by W'^ the equation becomes 

(52) 

This last identity is a useful reformulation of the problem for < r < 1, since in this 
case the two integral terms in the equation are positive. Then, eliminating both of 
them and integrating on [0,y], we get for each y >0 



W{y) - a 



T 



^^y^ - +7(erv/4_i) ■ (53) 
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For r > 1 it is more convenient to integrate by parts the first integral term in (I52p to 
obtain 



y ^' ^ j_- T y 



W 



i-r 1 ry 



w{z)dz + je-^y' 



+ 



T 1 

WW^ Jo 



W\z)Az = 0. 



Again, eliminating the two positive integral terms and multiplying by e 
obtain 

(e(^-i)./4^)'>le-/4. 
After an integration on [0, y], this gives 

(r-l) y/4^ > i + 1 -^-yl^ 



(r-l)a/4 



S{y)< 



ay 



ey/4 + a(e!//4-l) ' 
Summarizing, estimates (|53p and (|54|l read 

S{y) < ay g{y;a,T) , 

for all J/ > 0, a > 0, r > 0, where 

if < r < 1 , 



(54) 



(55) 



g{y-a,T) ■- < 



(r + a) e^y/'^-a 

1 

, (l + a)e2'/4-a 



(56) 



if r > 1 . 



As an important consequence of (|55[) - (|56|) . for any r > 0, 5 is bounded uniformly with 
respect to a > 0: 

min{r, l}y 



Siy)< 



^{T.l}y/i „ I 



(57) 



for each y > 0. Such an estimate does not follow from (|45[) - (|46|l . 

Estimate (|55|l is better than estimate (|45|l for r = 1. For r 7^ 1, this depends on the 
values of r and a. Therefore, it is interesting to reproduce the computations giving (|51|l 
by using the function g instead of h. For t > 1 and each j/ > 0, we obtain 



y 4 

g{z; a, r)dz = — log 
a 



[(l + «) 



and from equation (148 



e2'/4 

<t>{y )>a 

[{l + a)cy/^-a] 



This gives, for a > and r > 1, 



M(a,r) ^ 4 a 



27r 



a+ 1 



(58) 



(59) 
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Such a lower bound is definitely worse than (|5ip for large values of t or, to be precise, 
as soon as /(r) < log(a + l)/(2a). On the other hand, for r < 1, we have 

f g{z;a,T)dz = -log [(- + l) e'^^/* - -] --y, 

and again from equation (|48p , 

0'(y) > ae-^/S"* ^ . 

Finally, it holds true that, for a > and r < 1, 

Mhiilyl^, (60) 

To conclude, integrating (|58|l on [0, y] and using estimate (|55p gives us, for any r > 1, 

^'^^ - (1 + i) ef/4 _ 1 - (1 + i) ey/4 _ 1 - ^''^> 

for each y > 0, which is a good approximation of S and near the origin since it takes 
into account the condition ^'(O) — (/>' (0) = a. Moreover, (|49p is improved and replaced 
with 

^ M(a,T) (, _-y/i\ 4(c 



for any r > 1 and i/ > 0. 

Remark 5 As an additional consequence of the above estimates, we observe that 

a r g[z- a, r) d. = 4 log (l + .^e" '"'"{l-} ) 

J 

converges as y — > oo , so that 



exp 



g{z\a,T)Az 



2 

= (1 + 



min{l,r} 

According to (1411) . (|48l) and (|55p . we find the estimate 

§=/(a,r)= lim a exp ( -- dz | > a f 1 + . f, -. Y^, 

which, taking into account the change of parametrization s = log(2a), refines the 
estimate lim^^+oo c"(s) = found in 15 and our estimate (|25p . Here we use the 
notations of Remark [T] 
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4.3 New upper bounds 



Using the previous estimates on S and an argument in 12 , we can improve on the upper 
bound in (|43)) . Let 

oo if < T < i , 



if ^ <T<1, 



if r > 1 . 



Proposition 1 For any r > 0, if a < max{j(r), 1}, then M{a,T) < 87rmin{l,a}. 

The above estimate gives us a nonoptimal set of parameters (a, r) that guarantees 
M{a,T) < Stt. It is interesting to notice that \\xa^_^fj^2)+ Ji''') = 

Proof Let M = M(a,T). From the identity 

~^{^) - {^'l>iy)^'(y) + ^y^"iy)) ^y, 

and 4,y<f)" = -~y(f) — 2(j}S which follows from (|29p . we have, after an integration by 
parts and using (I30|l . 




Hence we have ^ < 4 if 

^ < 1 (61) 

for each y > Q. From (|57|) it follows that S{y) < 4 for all y > 0, for any r > and 
a > 0: the above sufficient condition (|61|l is satisfied whenever r < 1/2. For r > 1/2 
we have to use one of the previous upper bounds for 5*. 

(a) Using (|45p . we have 

l-ls{y)>l-2a-^{e-y/'-e--y/') 

for any r 7^ 1 and each y > 0, and condition (jGlf) is satisfied if 

^ . 1 r- 1 1 1 ^ 

a < mm 7- 7- = - r ^-1 . 

~ y>0 2 r e^y/"^ - c^'^ y/'^ 2 

For T — 1, using (|45p as before (or by continuity of the previous argument as r ^ 1), 
we similarly obtain 

^ . 2e^/4 e 
a < mm = - . 

y>0 y 2 

(b) Using (|55p . we have for r > 1 and each y > 
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Then condition (|6ip is satisfied if 

2e^/4 el- 

a < mm 



0<y<y ry-2(e3^/4-l) 2r 



e 



1- 



where we take into account that ry ~ 2 — 1) < for y > y, y being the unique 

solution of the equation -^y+l^e^/^. Similarly, for ^ < t < 1 and each j/ > 0, we 
get 

7^ S{jj)-l< ^ , . 1. 

Then condition (|6ip is satisfied if 

2Te4'' e 2T 

a < mm 



0<ry<y r^y -2[eiy ~1) 2t - e^'^ 

Comparing the results obtained in (a) and (fe), the proof of Proposition [T]is completed. 

□ 



5 Proofs 

This section is devoted to the proof of Theorems [51 13] and 3] As a byproduct of these 
results, we obtain Theorem [1] 



5.1 Proof of Theorem E] 

Given any fixed (a, r) £ , the local existence issue of the (singular) system (|29|l - (|30|) 
with initial conditions can be solved using a fixed point argument applied to the 
operator 

T[<5](j/) =ae-^/*- ie-^/'' T i e^^-^) <^(z) f e^l^<l>i£) d^, 
2 Jo ^ Jo 

defined on the complete metric space Xa := £ C[0, 2/a] : ^(0) = a, < ^(y) < a, 
< J/ < J/a} endowed with the usual supremum norm. Indeed, an appropriate choice 
of j/a gives that T maps Xa into Xa and that T is a contraction. If = ^, it 

is then enough to define (^{y) := ${z) Az and 5'(y) ~ (.-'^vl'^ <l>{z) dz in 

order that (0,5*) is a solution of JSHI-O, 132} with G C^[0,ya] n C^(0,ya] and 
5 G C^fO, j/a]. The continuation of the local solution to a global one is standard since 
system (|29|) - p0[) is no more singular away from the origin and solutions are locally 
bounded on R+ by the estimates of Section [4] 

The fact that cj> G C^[0,oo) follows from (gH)) and limy^o+ S{y)/y = a = S'(0). 
Estimates (|36} have been proved in Section [4] 

Finally, uniqueness of global solutions of (|29p - p0|l . (I33|l is a consequence of the 
contraction property of T and the Cauchy-Lipschitz theorem. 

Concerning the continuity of the map a G R+ {(j},S), let us denote by {(pi, Si) 
the solution associated to the shooting parameter Oj, i = 1,2. Following T2] we have 

|log,^i(y)-log02(y)l < |logai-loga2l + i T ^ - 52(^)1 d^ (62) 

^ Jo ^ 
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and 



\Si{y) - S2{y)\ < e--y/* r e(--l)-/4 I e'/U'liz) - e'^U2{z) I dz 

^gmax{logai.loga2}g-ry/4 f\{r~l)./A^ log (2) - log 0^ (2) | d^ , (63) 

^0 

where the decreasing monotonicity property of the function e^^* has been used 

in the last inequality. Plugging (|63} into ^ and denoting C = ^^^^{yog,ai,loga2} ^ 
obtain 



I log 01 (y) - iog(?!>2(y)| 



< |logai-loga2i + | r^e-^'^' / e^^-^) log 0i (C) - log </>2(C)l dC dz 







< j log ai - log a2 1 + I ^ 1 log 4>'i (C) - log <^2 (C) 1 /(C) dC , (64) 

where /(C) = e^^-^) /^°° i e"^ dz. Next, / G L^{Q,oo) with /o^/(C)dC = 
4 ~ 4/(r). Therefore, the Gronwall lemma applied to (|64|l gives us 

I log 0i(y)- log 02(2/) I < |logai-loga2|e^^o'/(C)dC < j log ai -log a2 I "^^^^ (65) 

Estimate (|65|) implies the continuity of the map a t-^ (j>' . The continuity of the maps 
ai-^ S and a 1-^ cj> follows by ((63l)-(l65j and by the identity (p{y) = S{y) + J S{z) dz 
respectively. Finally, the continuity of a M follows by ^ — J S{y) dy; see 
Section 0] for more details. □ 



5.2 Proof of Theorem |3] 

The existence of a sequence of positive self-similar solutions {u^.,Vf^) corresponding to 
a positive sequence {a^,} is an immediate consequence of the existence of a positive 
solution {4>k,Sk) of (f29)) - (f30)) . (f33)) by Theorem [2l Indeed, it is sufficient to define 

«fc(r) = 24(r2) and ^,(r) = i ^ dz , 

as follows from (|26|l and (|35p . Moreover, G C^[0, 00) and v^. £ C^[0,oo). Whenever 
afc — » 00, the limit ||wfc||ioo(K2-) 00 follows from ||'ffc||i:,oo(R2) = i'fe(O) and H24p . 

Next, let us define Mj, := ||wfe|lii(R2). From the estimates of Section |4j the sequence 
{Mj;} is bounded from above (by a constant depending on r), and there exist two 
subsequences, still denoted M/^ and Ufc, such that Mj, a and uj, ^ uSq. The delta 
measure is centered at ^ = since mj;(0) — 2aj^. Actually a — Sn for any t > 0, as an 
immediate consequence of the identity obtained in the proof of Proposition [1] 

Hence we have 

(§)'-"(§) -;/.."'«'(*^Kft-*"«i')-iifi^) <«) 

Letting k ^ 00 and observing that: 
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(i) — SkY ~ T Sk/4: is bounded by (|57p . uniformly with respect to aj. ^ oo, 

(ii) ui^{r) — 2 0^(r^) is uniformly decaying (with respect to aj. — > oo) for large values 
of r and limr— >oo sup^ /|^|>r"'=(^) 1^1^ '^^ = 0, as a consequence of (|40)) and (f57)) . 

we obtain that the right hand side in (I66|) converges to 0. On the other hand, a is 
necessarily positive by (|59|l and (|60|l . which proves that a = Stt. □ 

Remark 6 Let us observe that the identity 

4 M + 2 / u(0 V«(C) iAi- [ u{£) d^ = 

follows from equation Q multiplied by |^|^ and from the integrability of u given by 
([6]) and ([8]). Mimicking a standard computation for the parabolic-elliptic Keller-Segel 
system by writing v = —-^ log(-) *u~\- v, the above identity reads 

4M-^ + 2f m(0 VC(0 ■ C - / |Cl\(C)de = 0. 

2 71" JR2 JR2 

See for instance [5l[8] for more details. Therefore, we have found that Vw(5) • ^ = 
"/■(ICI^) ~ 'S'(I^I^) > 0. This is consistent with the fact that, from equation (|lip . one 
easily finds that (l){r'^) - S{r^) = v'{s) ds. 

5.3 Proof of Theorem H 

For any fixed r, let us define M* (r) = sup^^^Q M{a, r). Since M is bounded from above 
with respect to r, uniformly in a, continuous with respect to a, such that M(0, t) — 
and lima— >oo M(a, r) Stt, M* (t) is well defined and finite. The theorem is then a 
straightforward consequence of Theorem [2] and Proposition [1] □ 

6 Numerical results 

In this section, we numerically illustrate the above results. In particular, we show the 
existence of positive forward self-similar solutions with mass above 8 n and their multi- 
plicity when r is large enough. We follow two diff'erent approaches: first the formulation 
(|15p - (|16p . and then the cumulated densities formulation based on (|29p - (|30p . 

6.1 Bifurcation diagrams 

The computations giving rise to Figs. 1 and 2 are based on the parametrization provided 
by (|15p - (|16p . Numerically, one has to be careful with the origin and solve psp on the 
interval (e, oo) with the initial conditions 

w(e ; s) = s — — e* and w' {e ; s) = —- e e'^ , 

obtained by the Taylor expansion at e > 0, small enough, thus dropping higher order 
terms in e. Observe that by (|15p w"{0;s) — —e^/2. In case of Fig. 2, one has to 
compute M(s), which is given by ([T7|l . by solving M'(r) = 27re"'^'''*^ e"'' r with the 
approximate initial condition M(e) = yre^e*. 
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In Fig. 2, we recover that M{s) ^ Stt as s ^ oo. Moreover, for r large enough, 
there are two solutions corresponding to a given M larger than 8n, with M — Sir not 
too large. Since it is of interest to decide for which values of r solutions may have mass 
larger than Svr, the small rectangle in Fig. 2 (left) is enlarged in Fig. 2 (right). 

It can be numerically checked that solving the equations on (e, rmax) with rmax = 10 
gives a good approximation of the solution. Furthermore, here we took e = 10~^ and 

s e [-10,20]. 



2 5 ~'-2°o°°°^^ 












/ 






'5 



Fig. 1 The set of all positive solutions of Ava + f S ■ Vt)^ + ac^'e'^^^^^^ = in Cg(M^), 
where a = <t(s) = e"''-°°''*\ is represented by the multivalued diagram s i— > (log cr, log iio- (0)) 
for T = 10", a = —2, —1, . . . , 3. Recall that the solutions v^^ are radial and decreasing so that 
Va{0) = ll'^CT llioo (•]j2) . We observe that maxlog a"(s) appears as an increasing function of r. 




Fig. 2 Le/t.- The set of all positive solutions of Z\»;CT + §C-VDCT + o-e"<'e-l^l^/'' = 0inC^(R2)is 
now represented by the diagram s i— > (log(l+M(s)), logtio-(O)) for t = 10", a = —2, —1, . . . , 3. 
We observe that max M(s) appears as an increasing function of r. 

Right: The plot is an enlargement of the rectangle of Fig. 2 (left), with r = 0.60, 0.62, 0.64, 
. . . , 0.90. Numerically, the first solution with mass larger than Stt appears for r G (0.62, 0.64), 
which is far below the bound found in Section [2] This is not easy to read on the above figure, 
but it can be shown graphically by enlarging it enough. 
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6.2 Cumulated densities 

Plots and bifurcation diagrams of forward self-similar solutions can be computed in the 
framework of cumulated densities (|29p - p0p . (|33p . However, again one has to be careful 
with the singularity at the origin. As above, since for e > small enough, 5' ~ 0' ~ a 
on (0, e) and so 

S{y) = ay + 0{e') and ^"(j/) ~ (1 + 2a) + 0(e) , 

we practically solve (|29p - p0|l on (e,ymax) with the initial data 

</)'(e) = a — ^ (1 + 2 a) e , (?!)(e) = ae — ^ (1 + 2 a) and S{e)=ae. 
4 8 

for any y £ (0, e). Obviously, having fixed e > 0, one has to take a in such a way that 
(/>'(e) — a = o(a). Here, we choose e = 10^^. Finally, we shall approximate M from 
below by 0(j/max) with i/max large enough. Figs. 3 and 4 correspond to the cases r = 0.1 
and T — 10 respectively. For r — 0.1, the value Stt for the total mass is achieved only 
asymptotically in the limit a ^ oo. For t — 10, self-similar solutions with mass M 
larger than 8 vr exist for a large enough. Finally, Figs. 5 and 6 show the total mass as 
a function of a and t. 
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Fig. 5 Left: The value of mass 4>{oo) = M(a, r) / (2 n) in the logarithmic scale as a function of 
a, for T = 0.1 fc-^ with fc = 1, 2, . . . , 10. Right: An enlargement around the value M{a, t)/(2 tt) = 
4 in the logarithmic scale as a function of a, for t = 0.50, 0.55, 0.60, . . . , 1.00. 
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Fig. 6 The value of the maximal (in terms of a) mass 0(oo) = M»(T)/(2-7r) as a function 
of T. Numerically, the first solution with mass larger than Stt appears for t S (0.62,0.64), as 
already noticed at the level of Fig. 2 (right). This is again not easy to read on the above figure, 
but it can be shown graphically by enlarging it enough. 



7 Conclusions 

Self-similar solutions are much more than an example of a family of solutions. The 
experience of various nonlinear diffusion equations shows that they are hkely to be 
attracting a whole class of solutions, although this is still an open question for the 
parabolic-parabolic Keller-Segel model with large mass (see [14] for a result for small 
mass solutions). It is quite reasonable to expect that well chosen perturbations of these 
solutions asymptotically converge in self-similar variables to the stationary solutions 
we have found. This actually raises a much more interesting question, which is how 
to determine the basin of attraction of these self-similar solutions and to understand 
where is the threshold between solutions for which diffusion predominates and solutions 
which aggregate. Clearly, it is not going to be as simple as in the parabolic-elliptic case, 
where a single parameter, the total mass, determines the asymptotic regime. We can 
conjecture that blowup occurs for mass large enough and even, maybe, as soon as the 
total mass of the system is above 8 tt if initial data are sufficiently concentrated. 

The model considered in this paper is by many aspects ridiculously simple. See, 
for instance, [9] to get a taste of the variety of the nonlinearities that make sense 
even for a rather crude modelling purpose. Still, these models, in limiting regimes, 
asymptotically exhibit scaling properties similar to the ones of the parabolic-parabolic 
Keller-Segel model considered here. Therefore, we believe that the information gath- 
ered above, together with the methods that have been introduced, for instance, the 
cumulated densities reformulation of the model, should definitely be some valuable 
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piece of information in the study of the asymptotic behaviors of the equations used in 
chemotaxis. 
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